Practality in complex networks: critical and supercritical skeletons 
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Fractal scaling — a power-law behavior of the number of boxes needed to tile a given network with 
respect to the lateral size of the box — is studied. We introduce a new box-covering algorithm that 
is a modified version of the original algorithm introduced by Song et al. [Nature (London) 433, 392 
(2005)]; this algorithm enables effective computation and easy implementation. Fractal networks 
are viewed as comprising a skeleton and shortcuts. The skeleton, embedded underneath the original 
network, is a special type of spanning tree based on the edge betweenness centrality; it provides 
a scaffold for the fractality of the network. When the skeleton is regarded as a branching tree, it 
exhibits a plateau in the mean branching number as a function of the distance from a root. For non- 
fractal networks, on the other hand, the mean branching number decays to zero without forming a 
plateau. Based on these observations, we construct a fractal network model by combining a random 
branching tree and local shortcuts. The scaffold branching tree can be either critical or supercritical, 
depending on the small-worldness of a given network. For the network constructed from the critical 
(supercritical) branching tree, the average number of vertices within a given box grows with the 
lateral size of the box according to a power-law (an exponential) form in the cluster-growing method. 
The critical and supercritical skeletons are observed in protein interaction networks and the world- 
wide web, respectively. The distribution of box masses, i.e., the number of vertices within each 
box, follows a power law P m (M) ~ A/ - ' 7 . The exponent rj depends on the box lateral size Ib- For 
small values of is, rj is equal to the degree exponent 7 of a given scale-free network, whereas n 
approaches the exponent r = 7/(7 — 1) as Ib increases, which is the exponent of the cluster-size 
distribution of the random branching tree. Finally, we study the perimeter H a of a given box a, 
i.e., the number of edges connected to different boxes from a given box a as a function of the box 
mass Mb, a- It is obtained that the average perimeter over the boxes with box mass Mb is likely to 
scale as (H(Mb)) ~ Mb, irrespective of the box size Ib- 
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I. INTRODUCTION 

Fractal scaling recently observed [| in real-world scale- 
free (SF) networks such as the world-wide web (WWW) 
0, metabolic network of Escherichia coli and other 
microorganisms |3j , and protein interaction network of 
Homo sapiens [4j has opened a new perspective in the 
study of networks. SF networks || are those that ex- 
hibit a power-law degree distribution Pd(k) ~ kr 1 . De- 
gree k is the number of edges connected to a given vertex. 
Fractal scaling implies a power-law relationship between 
the minimum number of boxes Nb((b) needed to tile the 
entire network and the lateral size of the boxes Ib, i.e., 



(1) 



where ds is the fractal dimension This power- 

law scaling implies that the average number of vertices 
(Mb ((b)) within a box of lateral box size £b scales ac- 
cording to a power law as 



(M B (l B ))~£i 



(2) 



Here, the relation of system size N ~ Nb{£b)(Mb(£b)) 
is used. This counting method is called the box-covering 



method. At a glance, the power-law fractal scaling (QJ is 
not consistent with the notion of small-worldness (SW) 
of SF networks. SW implies that the average number of 
vertices within a distance £c from a vertex scales as 



(M c (£c)) ~ e tc ' to 



(3) 



where £q is a constant. This counting method is called 
the cluster-growing method. Here, subscripts B and C 
represent the box-covering and cluster-growing methods, 
respectively. The number of vertices M within a box is 
referred to as the box mass. This contradiction can be 
resolved by the fact that a vertex is (can be) counted 
only once (more than once) in the box-covering method 
(in the cluster-growing method) . 

Recently, it was suggested that the fractal scaling orig- 
inates from the disassortative correlation between two 
neighboring degrees Q or the repulsion between hubs 
H . Moreover, we showed that the fractal network con- 
tains the fractal skeleton [T(| underneath it; this skeleton 
is a special type of spanning tree, formed by edges with 
the highest betweenness centralities 0, 0] or loads . 
The remaining edges in the network are referred to as 
shortcuts that contribute to loop formation. The skclc- 
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ton of an SF network also follows a power-law degree 
distribution, where its degree exponent can differ slightly 
from that of the original network ^j- For fractal net- 
works that follow fractal scaling Q, we have shown j| 
that each of their skeletons exhibits fractal scaling simi- 
lar to that of the original network. The number of boxes 
needed to cover the original network is almost identical to 
that needed to cover the skeleton. Thus, since the skele- 
ton is a simple tree structure, it is more useful than the 
original network for studying the origin of the fractality. 

It was shown [| that the skeleton of the fractal net- 
work exhibits a non-dying branching structure, referred 
to as a persistent branching structure hereafter. A skele- 
ton can be considered as a tree generated in a branching 
process |l4| starting from the root vertex. This mapping 
can be applied to any tree. If a branching process occurs 
in an uncorrelated manner, the branching tree obtained 
from it exhibits a plateau, albeit fluctuating, in the mean 
branching number function Ti(d), which is defined as the 
average number of offsprings created by vertices at a dis- 
tance d from the root. Actually, the plateau is formed 
when h(d) is independent of d; this is denoted as n for fu- 
ture discussions. The branching tree structure obtained 
from the random branching process is known to be a frac- 
tal for the critical case 0,[H|, where the mean branching 
rate is (n) = 1. Here, (n) is defined as 



(n) = 2J nh r, 



(4) 



n=0 



where b n is the probability that a vertex will produce n 
offsprings in each step. Thus, (n) = n for the random 
branching tree. The fractal dimension of the SF branch- 
ing tree generated with branching probability b n ~ n _7 
is given by [Til Hr| 



dn = 




2) for 2 < 7 < 3, 
for 7 > 3. 



(5) 



Thus, the presence of the skeleton as the critical branch- 
ing tree served as a scaffold for the fractality of the fractal 
network. Then, the original fractal network is a dressed 
structure to the skeleton with local shortcuts; the number 
of shortcuts is kept minimal in order to ensure fractality. 
This idea is demonstrated by observing that the number 
of boxes in the fractal scaling for an original frac- 
tal network is similar to that of its skeleton When 
(n) > 1, we will show that a supercritical branching tree 
is also a fractal from the perspective of the fractal scal- 
ing Q j although exponential relation in the average box 
mass holds in the cluster-growing method. The su- 
percritical branching tree also exhibits a plateau in n(d). 
Based on these observations, we define a network to be a 
fractal if (i) it exhibits a power-law scaling Eq. (1) in the 
box-covering method and (ii) its skeleton is also a fractal 
with the persistent branching structure, i.e., a plateau 
exists in the mean branching number function n{d). 

Based on these findings, we introduced a fractal net- 
work model by incorporating the random critical branch- 



ing tree and local shortcuts |9(. In this paper, we will 
show that the model can also be generalized for the super- 
critical branching tree, thereby facilitating a better un- 
derstanding of fractal networks. For example, the model 
based on the supercritical branching tree can explain the 
subtle coexistence of SW and fractality as observed in 
the WWW. 

Due to the heterogeneity of the degree distribution in 
the SF networks, the distribution of box masses is also 
nontrivial. It exhibits a power-law tail with exponent rj, 



Pm{M B ) 



(6) 



in the box-covering method, whereas it exhibits a peak 
at a characteristic mass in the cluster-growing method. 
In this paper, we perform a detailed analysis of the 
real-world fractal network as well as the fractal network 
model, thereby showing that the box-mass distribution 
in the box-covering method for the fractal network can 
be explained by the branching dynamics. The exponent 
77 of the box-mass distribution is related to the exponent 
t that describes the size distribution of random branch- 
ing trees 01 ■ I n particular, for the critical SF branching 
tree, r is known to be []J| 



j7/(7-l) 
T =\3/2 



for 2 < 7 < 3, 
for 7 > 3. 



(7) 



The same value of r can be derived for the supercritical 
SF branching tree; however, the power-law scaling be- 
havior is limited to a finite characteristic size depending 
on (n) and 7 0. Thus, the cluster-size distribution fol- 
lows a power law for both the critical branching tree and 
the supercritical branching tree up to the characteristic 
size. 

In the first part of this paper, we present the frac- 
tal property of real-world networks and the model of the 
fractal network in detail as well as a further analysis of 
our previous work 0. Initially, in Sec. II we introduce 
a modified version of the box-covering method employed 
in this paper, following which we present the fractal scal- 
ing (Sec. Ill) and the mean branching number analysis 
(Sec. IV) for an extended list of complex networks, real- 
world and model networks. In the later part, we provide 
a general description of a model of fractal SF networks 
including the supercritical branching tree and study its 
property in detail in Sec. V. In Sec. VI, we examine the 
average box mass and the box-mass distribution for the 
fractal networks in the box-covering and cluster-growing 
methods. The average perimeter of a box as a function of 
the box mass is studied in Sec. VII. The summary follows 
in Sec. VIII. 



II. THE BOX-COVERING METHOD 

The fundamental relation of fractal scaling is based 
on the procedure referred to as the box-covering method 
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FIG. 1: (Color online) Schematic illustration of the box- 
covering algorithm introduced in this work. Vertices are se- 
lected randomly, for example, from vertex 1 to 4 successively. 
Vertices within distance £b = 1 from vertex 1 are assigned to 
a box represented by the solid (red) circle. Vertices from ver- 
tex 2, not yet assigned to their respective box are represented 
by the dash-dot-dot (black) closed curve, vertices from vertex 
3 are represented by dash-dot (green) circle and vertices from 
vertex 4 are represented by the dashed (blue) ellipse. 

PJ that calculates the number of boxes Nb needed to 
cover the entire network with boxes of lateral size Ib- 
This is analogous to the box-counting method normally 
used in fractal geometry . Song et al. introduced a 
new definition of the box applicable to complex networks 
such that the maximum separation between any pair of 
vertices within each box is less than £g. However, this 
particular definition has proved to be inessential for frac- 
tal scaling. Throughout this study, we utilize a different 
version of the box-covering method introduced here; this 
method involves sequential steps of box covering, thereby 
providing an easy implementation: 

(i) Select a vertex randomly at each step; this vertex 
serves as a seed. 

(ii) Search the network by distance Ib from the seed 
and assign newly burned vertices to the new box. 
If no new vertex is found, do nothing. 

(iii) Repeat (i) and (ii) until all vertices are assigned to 
their respective boxes. 

The above method is schematically illustrated in Fig. ^ 
It should be noted that vertices can be disconnected 
within a box, but connected through a vertex (or ver- 
tices) in a different box (or boxes) as in the case of box 2 
as shown in Fig.^ On the other hand, if we construct a 
box with only connected vertices, the power-law behav- 
ior Eq. is not observed. The box size £b used here is 
related to £s approximately as £s ~ 2(b + 1- A differ- 
ent Monte Carlo realization of this procedure ((i)-(iii)) 
yields a different number of boxes for covering the net- 
work. In this study, for simplicity, we choose the smallest 
number of boxes among all the trials. Although, this al- 
gorithm provides equivalent fractal dimension ds to the 
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FIG. 2: (Color online) Comparison of the box-covering meth- 
ods introduced by Song et al. Q (□) and in this paper (o). 
The results obtained from the two box-covering methods ap- 
plied to the WWW are plotted here. The two methods yield 
the same fractal dimension da ~ 4.1. 

one introduced by Song et al. [ij, it is easier to im- 
plement and is effective in computation. In Fig. wc 
compare the two box-covering methods applied to the 
WWW, demonstrating that the same fractal dimension 
is obtained. It should be noted that the box num- 
ber Nb we employ is not the minimum number among 
all the possible tiling configurations. Finding the actual 
minimum number over all configurations is a challenging 
task. However, in this paper, we focus on the problem 
of the fractal scaling within the framework of the box- 
covering algorithm introduced above. 

III. FRACTAL SCALING ANALYSIS 

We present fractal scaling analysis for real-world net- 
works, which are listed in the left column of Fig. rcffig3. 
We first examine fractal networks such as (a) the world- 
wide web (WWW) , (b) the metabolic network of E. coli, 

(c) protein interaction network (PIN) of H. sapiens, and 

(d) of S. cerevisiae fl9l |. Next, non-fractal networks such 
as (e) the actor network [2jj, (f) coauthorship network 
[211 . (g) Internet at the autonomou s sy stem (AS) level 
|22j . (h) Internet at the router level [23j, and (i) power- 
grid of the USA [24| are studied. The characteristics of 
these networks are listed in Table I. Note that in the 
previous study by Song et al., the protein interaction 
network of S. cerevisiae was classified as a non-fractal 
network. In this work, we use a different dataset °f 
high-confidence protein interactions, for which the PIN 
is a fractal network. 

For the fractal networks (a)-(d), the original network 
and its skeleton exhibit the same fractal scaling behav- 
ior, and the respective statistics of the numbers of boxes 
needed to cover them are almost identical as shown in 
the left column of Figs. 3(a)-(d). The fractal dimensions 
for these networks are measured to be « 4.1, 3.5, 2.3, and 
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TABLE I: Properties of real-world networks studied in this work. For each network, the number of vertices N, the average 
degree (k), the assortativity mixing index r, the average separation (d) of all pairs of vertices, and the maximum separation 
dmax among all pairs of vertices are tabulated. 



Name 


N 


(k) 


T 


Id) 


d 


category 


World-wide web 


325729 


6.7 


—0.05 


7.2 


46 


fractal and SW 


Metabolic network of E. coli 


2859 


4.8 


-0.16 


4.7 


18 


fractal and SW 


PIN of H. sapiens 


563 


3.1 


-0.14 


6.9 


21 


fractal but not SW 


PIN of S. cerevisiae 


741 


4.7 


0.41 


10.8 


27 


fractal but not SW 


Actor network 


374511 


80.2 


0.22 


3.7 




non-fractal and SW 


Coauthorship network (cond-mat) 


13861 


6.4 


0.16 


0.6 


18 


non-fractal and SW 


Internet at the AS level 


16644 


4.3 


-0.20 


3.7 


10 


non-fractal and SW 


Router network 


284805 


3.2 


-0.01 


8.8 


30 


non-fractal and SW 


Power grid of the USA 


4941 


4.9 


0.06 


8.5 


17 


undetermined and SW 



2.1 for (a) the WWW, (b) the metabolic network, PIN of 
(c) human, and that of (d) yeast, respectively. A power- 
law behavior is not observed for non-fractal networks, in 
which Nb{(-b) decays faster than any power law. We 
also study the fractal scaling for a randomly spanning 
tree of each network, which is constructed from edges 
that are randomly selected from the original network to 
form a tree. Since edges are selected randomly, the de- 
gree distribution of the ori gina l network is conserved in 
the random spanning tree (25J. The random spanning 
tree is fractal irrespective of the fractality of the ori gina l 
networks; this follows from the percolation theory |26j . 
Thus, the random spanning tree follows a power law in 
the fractal scaling. 

For non-fractal networks, the box number Nb does not 
follow a power law with respect to the box lateral size 
Hb ■ The statistics of box number of skeletons differ sig- 
nificantly from those of their original network; however, 
there are some exceptions in our examination. For the 
Internet at the AS level [Fig. EUg)] , n b for the original 
network and skeleton exhibits similar behavior; however, 
they do not follow a power law. On the other hand, in the 
power grid [Fig. Illlf il]. the fractal scalings of the origi- 
nal network and the skeleton exhibit the same power-law 
exponent; however, the box numbers for the original net- 
work and the skeleton differ significantly. Although the 
fractal scaling exhibits a power-law behavior, this fractal- 
ity is not obvious, because the network size is too small 
for checking if a plateau is intrinsically formed in n{d). 
Thus, the fractality cannot be classified. 

Next, we examine fractal scaling in the following net- 
work models (a) the Barabasi- Albert model with the de- 
gree of incident vertex m = 2 [jj, (b) Barabasi- Albert 
model tree with m = 1, (c) static model 0, (d) geomet- 
ric growth model [27]], and (e) deterministic hierarchical 
model in Fig. lIIII The network models considered do 
not obey the power-law fractal scaling; therefore, they are 
not fractals. 



IV. MEAN BRANCHING NUMBER ANALYSIS 

We present the mean branching number (MBN) anal- 
ysis for the skeleton and random spanning trees of each 
network considered. We define the MBN function n(d) 
as the mean number of offsprings of each vertex at dis- 
tance d from the root in a branching tree. For the fractal 
networks [Figs. 3(a')-(d')], both the skeleton and the ran- 
dom spanning tree exhibit a plateau in MBN, a signature 
of a persistent branching structure. For random spanning 
trees, the location of the plateau is distinctly obtained as 
1. With regard to the skeletons, while the plateaus 
in MBN of the protein interaction networks [Figs. 3(c')- 
(d')] appear to be located around n pa 1, they cannot 
be located clearly for the WWW and metabolic network 
due to large fluctuations [Figs. 3(a')-(b')]. Such fluctu- 
ations may originate due to various factors such as the 
finite-size effect and the artificial choice of the root of the 
branching tree. The dynamic origin of the formation of 
real-world networks may well be more complicated than 
the purely random branching dynamics: Thus, nontriv- 
ial correlations may exist. Although the location of the 
plateau in MBN cannot be clearly determined in some 
cases, its presence is a distinct feature of the fractal net- 
works and is absent in non- fractal networks. 

For non-fractal networks, the MBN of the skeleton de- 
cays to zero without forming a plateau [Figs. 3(e')-(h'))]. 
This is because the skeleton of each non-fractal SF net- 
work belongs to the class of "causal" trees [2|j, where 
vertices closer to the root are likely to have larger de- 
grees. In such structures, MBN decreases steadily with 
the distance from the root; therefore, a plateau cannot 
be formed. This absence of a plateau in MBN is also ob- 
served in the skeletons of the network models shown in 
Fig. lIIII Note that even for non- fractal networks, the ran- 
dom spanning trees exhibit plateaus at n ps 1, confirm- 
ing their fractality independent of the underlying original 
network structure. 

Although the fractal skeleton provides a scaffold for 
fractality in fractal networks, the manner in which the 
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FIG. 3: Fractal scaling analysis (left column) and mean branch- 
ing number (right column) of real-world networks including frac- 
tal (a)-(d) and non-fractal (e)-(i) networks. For each network, 
the original network (o), the skeleton (V), and a random span- 
ning tree ( ) are studied. In (a)-(d), the straight lines, drawn 
for guidance, have slopes of —4.1, —3.5, —2.3, and —2.1, respec- 
tively. In (e)-(h), the fit for the exponential function in the case 
of the original network and a power-law fit for the random span- 
ning tree are shown. In the right panels (a')-(i'), the horizontal 
line at 1 is drawn for reference. 



shortcuts are placed in the network is also important for 
preserving the fractality. With regard to this, the previ- 
ous result of the length distribution of shortcuts is 
important. It is known that two types of shortcut length 



distributions exist [10|. In the first type, the shortcut 
length distribution decays completely with respect to the 
shortcut length. In the other type, the shortcut length 
distribution exhibits a peak at a finite length, which is 
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FIG. 4: Fractal scaling analysis (a)-(e) and mean branching 
number (a')-(e') of the SF network models. The symbols used 
are similar to those in Fig. lIIII except in (d). The different sym- 
bols in (d) represents different parameters in the model. In (a) 
and (c), the exponential fit for the original network and a power- 
law fit for the random spanning tree are shown. For (b) and (d) , 
which are tree networks, the exponential fit is shown. In (e), 
the exponential fit to the original network is shown. Note that 
in (e), the box numbers for the original network, the skeleton, 
and the random spanning tree all overlap. In panels (a')-(e'), a 
horizontal line at 1 is drawn for reference. 



comparable to the average separation of all pairs of ver- 
tices on each skeleton; this is indicated by an arrow for 
each network in Fig. Thus, in the former case, the 
shortcuts connecting different branches of the skeleton 
are rare; however, their contribution is considerable in 
the latter case. In the latter case, the network will be 
globally interwoven and loses the fractality as in the case 
of random SF networks. Indeed, we have found that the 
fractal networks exhibit the former behavior [Figs.EJa) 
(d)], while the non-fractal networks exhibit the latter 
[Figs. OJe)-(j)]. This indicates that the shortcuts in the 
fractal networks are mainly local. 

V. FRACTAL NETWORK MODEL 

The observation of a plateau in the MBN for the skele- 
ton of the fractal networks prompted the construction of 
a fractal network model based on a random branching 
tree. We construct the model by reversing the steps, fol- 
lowed thus far to reveal the fractality. We first construct 
a branching tree in which the branching proceeds stochas- 
tically with a prescribed branching probability b n . We 
choose b n to follow a power law with respect to n in or- 
der to generate an SF network. Then, the branching tree 



is dressed with local shortcuts as well as global ones. The 
global connection is introduced to observe the crossover 
from fractal to non-fractal behavior. The frequency of 
global shortcuts is an important parameter of the model. 
More specifically, we consider the branching probability 
b n , i.e., the probability to generate n offsprings in each 
branching step, as 

b n = \n-i ( 7 >2) (8) 

for n > 1, and 

CO 

b = l-J2 b n (9) 

71=1 

for n = 0. Then, the resulting tree network is an SF 
network with the degree exponent 7. In order to generate 
a critical branching tree, the normalization constant Z 
is set to be Z = £(7 — 1). Where, £(x) is the Riemann 
zeta function, which follows from the criticality condition 
(n) = nb n — 1. We can also generate a supercritical 
branching tree by setting Z = £(7 — l)/( n ) with (n) > 1. 

After we generate a branching tree, we dress it with 
shortcuts by increasing the degree of each vertex by a 
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FIG. 5: Length distribution of shortcuts. The length of a 
shortcut is defined as the shortest distance along the skeleton 
between the two vertices connected by the shortcut. The 
arrows indicate the diameter of the skeleton of each network. 
Here, n s (d) and tit are the number of shortcuts with length 
d and the total number of shortcuts, respectively. 



factor p and establishing the available connections be- 
tween vertices. This can be achieved either in a local or 
global manner. An additional parameter q is introduced 
to describe the frequency of global shortcuts in the net- 
work. The creation rule of the fractal network model is 
described as follows: 



(i) We start with a seed vertex from which n off- 
springs are stochastically generated with probabil- 
ity b n (n = 0, N — 1). Each offspring then gener- 
ates n branches with probability b n . This process 
is repeated until we obtain a network of desired size 
N. If the growth of the tree stops before attaining 
size N, we restart the branching procedure. 

(ii) Degree fcj of each vertex i is increased by a factor p 
such that vertex i obtains additional pki stubs for 
forming edges. From these stubs, qpki stubs are 
assigned to global shortcuts, while the remaining 
(1 — q)pki stubs are assigned to local shortcuts. In 
order to establish local shortcuts, we search vertices 
from the root. A vertex i that has at least one stub 
for local shortcuts is selected. Then, its connection 
partner is scratched from the closest vertices from 
the vertex i to a vertex j, having available stubs for 
local shortcuts and not yet connected to i, to form 
an edge between i and j. This process is repeated 
until all the stubs for local shortcuts are linked. 

(iii) Next, we choose two vertices i and j randomly; each 
of which has at least one stub for global shortcuts. 
We then connect them to form an edge if they are 
not already connected. This process is repeated 
until all stubs for global shortcuts are linked. This 
step is similar to the process used in the configura- 
tion model |3lj . 

Network configurations obtained by the model with 
7 = 2.3 and N = 164 are shown in Fig. HjJ A critical 
branching tree with (n) = 1 is shown in Fig. Eta); this 
tree is dressed by local [Fig.EJb)] generated with p = 0.5 
and q = and both local and global shortcuts [Fig. Etc)] 
with p = 0.5 and q = 0.02. A supercritical branching tree 
with (n) — 2 is shown in Fig. EJ^d); this tree is dressed 
by both local and global shortcuts generated with pa- 
rameters p = 0.5 and q = 0.02 in Figs. Etc) and (f), 
respectively. 

We examine the fractal scaling in the network model 
and the MBN for its skeleton. In the case of a network 
generated from a critical branching tree (with 7 = 2.3 
and N « 3 x 10 5 ) and dressed only by local shortcuts 
(with p = 0.5 and q = 0), 76% of all edges of the 
original branching tree are maintained in the skeleton. 
The branching tree and the dressed network exhibit frac- 
tal scalings with the same fractal exponent ds ~ 3.2 
[Fig. Eta)]. This value appears to differ from the theoret- 
ical value w 4.3 estimated from the formula JSJl. However, 
we notice that the measured value of the degree exponent 
of the dressed network is rather close to 7 = 2.4, although 
the branching tree is generated with parameter 7 = 2.3. 
Thus, the expected value is ds = 3.5. Therefore, the 
numerical deviation can be explained. The MBN of the 
skeleton of the dressed network displays a plateau around 
1 [Fig. Ett>)]. Moreover, when we introduce 1% of global 
shortcuts (p — 0.5 and q — 0.01) to the critical branch- 
ing tree, the box number Nb(£b) decays faster than any 
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FIG. 6: (Color online) Snapshots of the fractal network models, (a) A critical branching tree with 7 = 2.3 and TV = 164 created 
in step (i). (b) A network created by adding local shortcuts (blue) following the rule (ii) to the branching tree in (a). Parameter 
p = 0.5 is used. This network is still fractal, (c) A network created by adding global shortcuts (red) following the rule (iii) 
to the network in (b). Parameter q = 0.02 is used. This network is no longer fractal, but small- world, (d) A supercritical 
branching tree with mean branching number (n) = 2, which is fractal as well as small-world, (e) A dressed network to network 
(d) by local shortcuts (blue) generated with p — 0.5. The network is fractal as well as small-world, (f) A network created by 
adding global shortcuts (red) to the network in (e). q — 0.02 is used. In (a)-(f), the colors of each vertex represent distinct 
generations from the root. 



power law for large values of Ib [Fig. 13a)]; thus fractality 
is lost. Accordingly, in this case, the MBN of the skeleton 
decays to zero without a plateau [Fig.0[b)]. The critical 
value q c above which the network becomes non-fractal 
depends on the degree exponent 7, system size N, and 
number of shortcuts p. A more detailed analysis on this 
crossover behavior will be presented elsewhere |32| . 

The same analysis is performed for the model based 
on the supercritical branching tree with (n) = 2 [Fig. [S] . 
This tree with 7 = 2.3 displays a power-law fractal scal- 
ing with fractal exponent ds ~ 4.2 [Fig. Eta)]; how- 
ever, its MBN fluctuates heavily on and off about the 
expected value n = 2, while exhibiting persistent branch- 
ing [Fig. IHtb)]. Such a highly fluctuating MBN is sim- 
ilar to that observed in the skeleton of the WWW or 
the metabolic network [Figs. lTTlT a'Wb')]. When dressed 
only by local shortcuts (p = 0.5 and q = 0), the 
dressed network still exhibits a power-law fractal scaling 
[Fig. Efa)]. The MBN of its skeleton still exhibit large 
fluctuations [Fig. Etb)]; however, its mean n decreases 
from 2. With 1% of global shortcuts (q = 0.01), the 



fractal scaling exhibits a power-law behavior but with an 
exponential cutoff [Fig. |SJc)]. Interestingly, the MBN of 
its skeleton displays a plateau located around 1 with re- 
duced fluctuations [Fig. EJd)]. When we further increase 
the number of global shortcuts to 2% (q = 0.02), the 
MBN of the skeleton decays without a plateau, and the 
network is no longer a fractal [Figs. [S^c) -(d)]. 



The fractal network model studied here is a generaliza- 
tion of the previous model [t| generated by including the 
supercritical branching case. Using this model, we can 
reproduce the highly fluctuating behavior in the MBN 
observed in the WWW as well as the fractal scaling. We 
will also show that this generalization of the supercrit- 
ical branching facilitates a better understanding of how 
the SW and fractal scaling coexist and do not contradict 
each other in such systems. In the following, we investi- 
gate the properties of the fractal network using the fractal 
network model as well as the real- world fractal networks. 
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2 4 8 16 32 64 5 10 15 20 25 30 35 40 45 50 

Ib Distance from the root 

FIG. 7: (Color online) Fractal scaling analysis (a) and mean 
branching number (b) for the fractal models based on the 
critical branching tree with only local shortcuts with p = 0.5 
and 9 = (□), and with 1% of global shortcuts with p — 
0.5 and q = 0.01 (A). The bare critical branching tree is 
represented by (O)- The solid line in (a) is guideline with a 
slope of -3.2. The measured degree exponent is 7 ~ 2.4 and 
system size is N = 3 x 10 5 . 
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FIG. 8: (Color online) Fractal scaling analysis (a) and (c), 
and mean branching number (b) and (d) for the fractal models 
generated from a supercritical branching tree with (n) = 2, 
dressed by shortcuts. The data are for the bare supercritical 
tree (o) and dressed networks with p = 0.5 and q = (□), 
p = 0.5 and q = 0.01 (A), and p = and q = 0.05 (o). 
The solid lines in (a) and (c) are guidelines with slopes of 
-4.2 each. The degree exponent is 7 = 2.3 and system size is 
N = 1 X 10 5 . 



VI. SMALL-WORLDNESS AND BOX MASS 
DISTRIBUTION 

In this section, we study the average box mass (M(£)) 
as a function of box size £. The box mass is measured 
using two methods, the cluster-growing method and the 



box-covering method In the cluster-growing method, 
the box mass is defined as the number of vertices at 
a distance not greater than £c from a given vertex. 
Note that in the cluster-growing method, a vertex can 
be counted by more than one box, whereas in the box- 
covering method, it is counted only once. The cluster- 
growing method provides information on the SW of the 
network. The average box mass for the critical branch- 
ing tree grows with distance £c, according to a power law 
with the exponent ds, 

(M c (£c))~t d c B , (10) 

where ds is defined in Eq. J^J; this implies that the 
critical branching tree is a fractal. For the supercriti- 
cal branching tree, the average mass grows exponentially 
with increase in distance £q, 

(M c (£c))~(n) e °. (11) 

This relation is equivalent to Eq. (j2J), thereby suggest- 
ing that the supercritical branching tree is a small-world 
network. On the other hand, the average box mass in 
the box-covering method is determined by the fractal- 
ity. Since both the critical and supercritical branching 
trees are fractals, their box mass increases according to a 
power law. Thus, the analysis of the average box mass in 
the two methods will provide an insight into the interplay 
between SW and fractality. 

We also consider the distribution of box masses P m (M) 
for each method. It is known Q] that in the cluster- 
growing method, the box mass distribution exhibits a 
peak at a characteristic mass, while it exhibits a fat tail 
without a peak in the box-covering method. However, 
the origin of the power-law behavior of the box mass 
distribution for the fractal networks has not been under- 
stood clearly. Here, we present a detailed analysis of the 
box mass distribution, showing that the exponent of the 
power-law behavior depends on the lateral size of the box. 
When the lateral size of the box is large, the exponent 
of the box-mass distribution can be understood from the 
perspective of branching dynamics. 

A. Real-world fractal networks 

We first examine the average box mass (Mc{£c)) m 
the cluster-growing method for the original network and 
the skeleton of each real-world fractal network. For the 
WWW, we find that both the original network and the 
skeleton exhibit an exponential increase in average box 
mass with distance £c [Fig-EJa)]. Thus, the WWW is a 
small-world network and the skeleton of the WWW is a 
supercritical branching tree. For the metabolic network, 
while the original network is small- world with Eq. Ullfl . 
its skeleton appears to follow a power law, Eq. 1|1(J|) 
[Fig. Etb)]. Thus, the skeleton of the metabolic network 
is better described by a critical branching tree, although 
the original network is small-world. The difference be- 
tween the metabolic network and its skeleton probably 
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FIG. 9: Average box mass {Mc(£c)) in the cluster-growing 
method divided by the total number of vertices N, as a func- 
tion of the distance ic- Plots of (a) and (b) are drawn on a 
semi-logarithmic scale and those of (c) and (d) on a double- 
logarithmic scale, respectively. Filled and open symbols rep- 
resent the original network and the skeleton of each network, 
respectively. The solid lines for reference in (c) and (d) have 
slopes of 1.9 and 2.3, respectively. 



originates from the presence of core subnetworks in the 
metabolic network, wherein the vertices are tightly inter- 
woven through multi ple p athwa y s but are simplified into 
a tree in the skeleton [3(|, 13a. l34j . On the other hand, for 
the protein interaction networks, both the original net- 
works and the skeletons behave according to a power-law 
form of Eq. i|10|) . Therefore, the protein interaction net- 
works are not likely to be small-world and their skeletons 
can be regarded as critical branching trees [Figs. EJc)- 
(d)]. 

Next, we study the average box mass {Mb{Ib)} in 



TABLE II: Behavior of the average box mass of the fractal 
networks and their skeletons in the cluster-growing and box- 
covering methods. 





Cluster-growing 


Box-covering 




method 


method 


World-wide web 


Exponential 


Power law 


World-wide web (skeleton) 


Exponential 


Power law 


Metabolic network 


Exponential 


Power law 


Metabolic network (skeleton) 


Power law 


Power law 


PIN of H. sapiens 


Power law 


Power law 


PIN of H. sapiens (skeleton) 


Power law 


Power law 


PIN of S. cerevisiae 


Power law 


Power law 


PIN of S. cerevisiae (skeleton) Power law 


Power law 
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FIG. 10: Average box mass {Mb(£b)) divided by the total 
number of vertices N, as a function of box size £b in the 
box-covering method. Filled and open symbols represent the 
original network and the skeleton of each network, respec- 
tively. The solid lines (drawn for reference) have slopes of 
4.1, 3.4, 2.0, and 1.8 for (a), (b), (c), and (d), respectively. 



the box-covering method. For the fractal networks and 
their skeletons, the average mass (Mb(£b)) increases ac- 
cording to a power law with respect to box size Ib of 
Eq. J2J) , regardless of whether it is critical or supercritical 
in the cluster-growing method as shown in Fig. 1101 The 
fractal dimensions measured using the formula © are 
d B = 4.1,3.4,2.0, and 1.8 for the WWW (a), metabolic 
network (b), protein interaction network of H. sapiens 
(c), and S. cerevisiae (d), respectively. These values are 
comparable to the ones obtained from the fractal scaling 
©, which are d B = 4.1, 3.5, 2.3, and 2.1 for (a), (b), 
(c), and (d), respectively. The results of the average box 
mass for the real-world fractal networks are summarized 
in Table II. Non-fractal networks exhibit the exponential 
relationship 



(M B (lB))~exp(£ B /lo) 



(12) 



with a constant £q. 

The different behaviors of the average mass in the two 
methods, the cluster-growing and box-covering methods, 
originates from whether overlap between the boxes is 
allowed. Thus, studying the extent of overlap of the 
boxes during the tiling can provide important informa- 
tion. In this regard, we measure the cumulative frac- 
tion F c (f) of vertices counted / times or more in the 
cluster-growing method for the WWW in Fig. ^2 The 
cumulative fraction F c (f) is likely to follow a power law 
for small /, thereby indicating that the overlaps occur in 
a non-negligible frequency even for a small distance £c- 
The associated exponent decreases with increase in box 
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FIG. 11: (Color online) Cumulative fraction F c (f) of the ver- 
tices counted / times in the cluster-growing algorithm. F c (f) 
follows a power law in the small / region, where the slope de- 
pends on box size lc- However, for large values of /, the data 
largely deviate from the value extrapolated from the power- 
law behavior. Data are presented for lc = 2 (•), lc = 3 (■), 
and e c =5( ). 

size £c as the chances of overlaps increase. However, for 
large values of /, the large fraction of vertices counted 
exceed the frequency extrapolated from the power-law 
behavior. As opposed to a bounded distribution such as 
a Poisson-type distribution, the broad distribution of / 
implies that a significant fraction of vertices are counted 
more than once in the cluster-growing method. Such 
multiple counting due to overlap is excluded in the box- 
covering method. Due to this exclusion effect, the mass 
of a box in the box-covering method is significantly lower 
than that in the cluster-growing method. 

We study the box-mass distributions in the two meth- 
ods. As shown in 0, for the WWW, the box-mass dis- 
tribution in the cluster-growing method exhibits a clear 
peak [Fig. I12f a)] ; on the other hand, in the box-covering 
method, it exhibits a fat tail, following an asymptotic 
power law, 

P m (M B ) ~ Mp (13) 

[Fig. I12f bt] . We find that the exponent r/ depends on 
the box size Ib- For small tg — 1 or 2, it is found that 
X] is equal to 7; however, as Ib increases, r\ approaches 
the exponent r of the cluster-size distribution 0. This 
can be understood as follows. For small values of £b, 
the branching has not proceeded sufficiently to exhibit 
asymptotic behavior; thus, the box mass will simply scale 
with the degree of the seed vertex, which is selected ran- 
domly, yielding r\ = 7. This is most evident for £ B = 1. 
On the other hand, as £b increases, the box grows and 
its size governs the scaling. The growth of the box can 
be approximated by the SF branching tree with the ex- 
ponent 7, the size distribution of which follows a power 
law with the exponent r given by Eq. (0 , yielding 77 = r 
for large Ib ■ The numerical estimates of r\ obtained from 
the WWW are in reasonable agreement with the predic- 
tion as r) = 7 ~ 2.3 for Ib = 2 and 77 w 1.8 for Ib = 5 
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FIG. 12: (Color online) Box mass distribution in the cluster- 
growing method (a) and box-covering method (b) for the 
WWW. Data in (a) are for lc = 2 and those in (b) are for 
Ib = 2 (•) and Ib = 5 (■). The solid lines are guidelines 
with slopes of -2.2 and -1.8, respectively. 




FIG. 13: Average box mass versus box lateral size in the 
cluster-growing method (a) and the box-covering method (b) 
for the fractal network model constructed from a critical 
branching tree, dressed by shortcuts with p = 0.5 and q = 
(•) and p = 0.5 and q = 0.01 (■). Degree exponent is 7 = 2.3 
and system size is N = 3 x 10 5 . Solid lines in (a) have slopes 
of 3.3 and 2.4, respectively, and the solid line in (b) has a 
slope of 3.3. 

[Fig.12]. 



B. Fractal network model 

Here, we perform a similar analysis of the average box 
mass and the box-mass distribution for the fractal net- 
work model introduced in Sec. V. We first consider a 
network model based on a critical branching tree with 
mean branching number (n) = 1. For simplicity, we fix 
the parameters to be 7 = 2.3, JV = 3x 10 5 , and p = 0.5, 
while varying parameter q. When q is sufficiently small, 
i.e., q < 0.001, the model network exhibits a power-law 
scaling both in the cluster-growing method Eq. (|10fl and 
the box-covering method Eq.(j2J. The network thus re- 
mains as a fractal. However, for larger values of q like 
0.01, the fractal scaling breaks down and the average box 
mass increases exponentially as Eq. <|11|) in both methods 
[Fig.^5], i.e., the network becomes small-world. 

The behavior of the average box mass of the fractal 
model network based on the supercritical tree is inter- 
esting. Once a supercritical tree is generated, the model 
network is dressed by local shortcuts with p = 0.5. Then 
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FIG. 14: (Color online) Average box mass as a function of box 
size in the cluster-growing method (a) and the box-covering 
method (b) for the model network constructed from a super- 
critical branching tree, dressed by shortcuts with p = 0.5 and 
q = (•) and p = 0.5 and q — 0.01 (■). Degree exponent is 
7 = 2.3, and system size N = 1 X 10°. The solid line in (b) 
has a slope of 4.0, which is drawn for guidance. 




FIG. 15: (Color online) Box-mass distribution in the cluster- 
growing method (a) and box-covering method (b) for the frac- 
tal model network grown from a critical branching tree with 
7 = 2.3 and dressed by shortcuts generated with p = 0.5 and 
q = 0.001. The data in (b) are for £ B = 2 (•) and l B = 5 (■). 
Their slopes are -2.3 and -1.8, respectively. The system size 
is N » 3 x 10 5 . 




FIG. 16: (Color online) (a) Average box mass versus box size 
in the cluster-growing (x) and box-covering (■) methods for 
a bare critical branching tree with 7 = 2.6. Solid lines, drawn 
for guidance, have slopes of 2.3 (x) and 2.6 (■), respectively, 
(b) Box-mass distribution for the bare tree of (a) . Solid guide- 
lines have slopes of -2.8 for £b — 2 (•), -2.6 for £b — 5 (■), 
and -1.6 for l B = 32 (a). 
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FIG. 17: (Color online) Box-mass distribution in the cluster- 
growing method (a) and box-covering method (b) for the frac- 
tal model network grown from a supercritical branching tree 
with (n) = 2 and 7 = 2.3, which is dressed by shortcuts 
with the parameters p — 0.5 and q = 0. The box size in (b) is 
Im = 2 (•) and £m = 5 (■). Solid lines in (b) have slopes -2.3 
and -1.8, drawn for guidance. The system size is N « 1 x 10 5 . 



it simultaneously exhibits both an exponential increase 
in box mass in the cluster-growing method and a power- 
law increase in the box-covering method, as observed in 
the WWW. This coexistence persists when we introduce 
global shortcuts up to q = 0.01 [Fig. El- If 

we further in- 
crease q, the average box mass increases exponentially in 
both methods as the network loses fractality. Thus, the 
model network with supercritical branching tree and an 
appropriate number of local shortcuts can reproduce the 
small- world property of the average box-mass function as 
well as the fractality. 

Next, we study the box-mass distribution for the frac- 
tal network model. We restart the analysis with the 
model network based on the critical branching tree with 
parameters 7 = 2.3, N = 3 x 10 5 , and p = 0.5. As 
with the real-world fractal networks like the WWW, the 
box-mass distribution for the model network exhibits a 
peak at a finite mass in the cluster-growing method. The 
box-mass distribution in the box-covering method fol- 
lows an asymptotic power law with exponent rj. As ob- 
served for the WWW, we observe that the exponent rj 



depends on the box size. For q = 0.001, it is found that 
r) w 2.3 for small £b = 2, and rj w 1.8 for large £b = 5 
[Fig. US]. The latter value 77 « 1.8 is in agreement with 
r = 7/(7— 1) ~ 1.8 from Eq. Q. Such ^-dependent be- 
havior of the box mass distribution can also be observed 
for another value of 7, for example, 7 = 2.6 [Fig. I16|. In 
such cases, the behavior 77 = r appears for large values 
of £bi for example, £b = 32 for 7 = 2.6. 

Next, when the model network is constructed based on 
a supercritical branching tree (with 7 = 2.3 and (n) = 2) 
and is dressed by shortcuts (with p = 0.5 and q = 0.001), 
rj w 2.3 is measured for £b = 2, however, 77 w 1.8 for 
£b = 5 [Fig. El- The obtained value r] « 1.8 is again 
in agreement with the expected value r « 1.8 for 7 = 2.3. 



VII. PERIMETER OF A BOX 

The boundary of a fractal object is an important phys- 
ical quantity and is considered to be another fractal ob- 
ject. For example, the area of the spin domain of the 
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FIG. 18: (Color online) Plot of the average perimeter 
(H(M B )) versus box mass Mb for the WWW (a) and the 
model network with the parameters of 7 = 2.3, p — 0.5, and 
q — 0.0 (b). Data are for box sizes Ib = 2 (•), £b = 3 (■), 
Ib = 5 (a). The solid line (drawn as reference) has a slope of 
1.0 for both (a) and (b). 



Ising model at critical temperature, which corresponds 
to magnetization, is a fractal object, and the interface 
length of the spin domain is another fractal object, cor- 
responding to the singular part of internal energy |35l | . 
Moreover, a percolation cluster is a fractal object and 
its outer boundary, referred to as "hull," is also a frac- 
tal 

We define the perimeter H a of a given box a as the 
number of edges connected on one end to the vertices 
within the box a and on the other end to the vertices in 
other boxes. The perimeter H a is examined as a function 
of the box mass Ms,a of the box a. Then, we can define 
the average perimeter (H(Mb)) over the boxes with box 
mass Mb- We find that the following power-law relation- 
ship exists, 



(H(M B )) ~ M 



(14) 



The new exponent dn (the hull exponent for the fractal 
network) is analogous to the one used in the percolation 
theory jME^. 

The power-law relation (|Tlfl is tested for the WWW 
and the network model generated with 7 = 2.3 and 
(n) = 1. It appears that (H(Mb)) depends on the box 
mass Mb linearly, i.e., dn/ds ~ 1, irrespective of £b 
for the WWW [Fig. EI a )]; however, it depends on £ B 
weakly for the fractal network model. For Ib = 2 and 3, 
dn/ds ~ 1; however, for Ib = 5, du/ds is likely to be 
marginally smaller than 1 [Fig. llSf b)]. The linear behav- 



ior implies da = ds, and is observed in the connections 
between the percolation clusters near the critical point 
for SF networks [H. 



VIII. CONCLUSIONS 

Recently, it was shown that some SF networks exhibit 
fractal scaling, N b (£b) ~ where N b (£b) is the 

number of boxes needed to tile the entire network with 
boxes of size Ib- In this paper, we have introduced 
a modified version of the box-covering method, which 
makes implementation easy. The origin of fractal scaling 
is understood from the perspective of criticality and 
supercriticality of the skeleton embedded underneath 
each fractal SF network. By performing the analysis of 
the average box mass as a function of the box size for 
the box-covering and cluster-growing methods, and the 
mean branching number as a function of the distance 
from the root, we found that the skeleton of the WWW 
is a supercritical branching tree, while the skeletons of 
other biological networks such as the metabolic network 
of E. coli and the protein interaction networks of H. 
sapiens and S. cerevisiae are critical branching trees. 
Based on this observation, we constructed the fractal 
network model. The box mass is heterogeneous and they 
exhibit a fat-tailed behavior, P m (M) ~ M~ v . We found 
that the exponent r\ depends on the lateral size ^ of the 
box. When £b is small, r\ is equal to the degree exponent 
7; on the other hand, as £b increases, rj approaches the 
exponent t = 7/(7 — 1) for the cluster-size distribution 
of the branching tree; this can be predicted from the 
skeleton. Finally, we studied the number of edges that 
interconnect a given box and other boxes, forming the 
perimeter of a box, as a function of the box mass. It 
appears that the perimeter depends on the box mass 
linearly, and the perimeter exponent is equal to the 
fractal dimension. 
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